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We consider the dynamics of the disordered trap model, which is known to be completely out-of- 
equilibrium and to present strong localization effects in its aging phase. We are interested into the 
influence of an external force, when it is applied from the very beginning at t — 0, or only after a 
waiting time t w . We obtain a "non-linear Fluctuation Theorem" for the corresponding one-time and 
two-time diffusion fronts in any given sample, that implies the following consequences : (i) for fixed 
, times, there exists a linear response regime, where the Fluctuation-Dissipation Relation or Einstein 

■ relation is valid even in the aging time sector, in contrast with other aging disordered systems; (ii) 

f"^ | for a fixed waiting time and fixed external field, the validity of the linear response regime is limited 

, in time by a characteristic time depending on the external force; (iii) in the non-linear response 

regime, there exists a very simple relation for the asymmetry in the position for the one-time and 
the two-times disorder averaged diffusion fronts in the presence of an external force, in contrast to 
other models of random walks in random media. The present non-linear Fluctuation Theorem is a 
consequence of a special dynamical symmetry of the trap model. 
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I. INTRODUCTION 



The field of non-equilibrium statistical mechanics is extremely broad, because the non-equilibrium character of 
the dynamics may have various origins. For instance, a system can be out-of-equilibrium because it is driven in a 
non-equilibrium steady state by external boundary conditions or driving forces [ij, |2J, y| or because it remains forever 
non-stationary and presents aging, such as in coarsening dynamics Q or in disordered models 0, 0, • 

Although these various classes of systems usually give rise to different questions and different approaches, it seems 
nevertheless that a central goal in all these non-equilibrium systems has been to study the possible generalizations of 
the Fluctuation-Dissipation Theorem, which is an essential result of Equilibrium statistical mechanics : it relates the 
' response to an external perturbation with the fluctuations existing in the absence of the perturbation. 

On one hand, for non-equilibrium steady states, various generalizations of the Fluctuation-Dissipation Theorem 
have been proposed @,EJE|, as well as the Gallavoti-Cohen Fluctuation Theorem for deterministic thermostated 
dynamics under chaoticity assumptions. This Fluctuation Theorem, which involves the time-reversed dynamics, has 
been then extended to various stochastic dynamics, such as Langevin dynamics |l2j . jump Markov processes and 
q | diffusion processes |l3l ll lj and allows to obtain results in the non-linear response regime. 

On the other hand, in the field of aging dynamics of disordered models, the proposed generalization of the 
Fluctuation-Dissipation Theorem concerns the linear response regime and consists in introducing a Fluctuation- 
• i— i , Dissipation ratio to describe the modified relation between the two-time response function and the two-time cor- 
relation 0, B • The introduction of this FDT violation ratio for the linear response regime has been also very 
useful to characterize the coarsening dynamics in non-disodered systems [l5L Il6l Il7j . Very recently, a non-linear 
Fluctuation Theorem has also been proposed for the aging dy namics of disordered spin models |Tsj. where however 
the time-reversal operation used in non- aging systems (lTLll2t fl3l fl^ | has to be completely redefined. 

Another important generalization of the usual Fluctuation-Dissipation Theorem deals with the the non-linear 
response of systems driven far away from thermal equilibrium, but where the initial state before the application of an 
external force is the thermal equilibrium [T^ . l2(i| . (Note that the initial state in the other approaches described above 
is usually not arbitrary but the Gibbs-Boltzmann distribution in |l2j or the stationary measure in jl3j|). Remarkably, 
the Jarzynski non-equilibrium work relation |Ts| has recently become a very useful tool to reconstruct the free-energy 
from non-equilibrium single-molecule pulling experiments |2 1| . 

This short and very incomplete presentation of recent studies of response and fluctuations in non-equilibrium 
statistical mechanics already shows that the subject is very rich and very diverse, and that some unifying framework 
is still lacking in comparison with equilibrium statistical mechanics. 

In this paper, we are interested into the response of trap models on a hypercubic lattice in finite dimension d 
to an external force field. Trap models provide a simple phenomenological mechanism for aging [22I l23l 124^. The 
one-dimensional version d = 1 is moreover interesting on its own since it appears in various physical applications 
concerning transport properties in disordered chains |25l [26] or in the dynamics of denaturation bubbles in random 
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It s ag ing and sub-aging properties in the absence of a bias have been much studied, either by 
29ll30j| and by physicists |3lll32[. Recently, the response of the one-dimensional trap model was 



2 



studied via scaling arguments and numerical simulations in |33| . where various regimes were found depending on the 
relative values of the two times considered (t w ,t) and the external applied field /, the main results being that the 
Fluctuation-Dissipation Relation (or Einstein relation) of the linear response regime is valid even in the aging sector, 
whereas the response always become non-linear at long times. The aim of this paper is to show that these response 
properties are actually two direct consequences of a non-linear Fluctuation Theorem, that we establish by using a 
special symmetry of the master equation of the trap model. Another important outcome of the non-linear Fluctuation 
Theorem is a prediction for the asymmetry induced by the external force on averaged diffusion fronts. We will always 
discuss the d = 1 case to simplify the presentation, but it is clear that the results can be straightforwardly generalized 
in arbitrary finite dimension d. However, we will not discuss mean-field trap models, where previous studies on the 
Fluctuation-Dissipation relation have shown that the results depend on the observable |34[ , and on the choice of the 
functional form of the rates [35l l3l|. For the trap models on lattice in arbitrary finite dimension that we consider, 
there are no such ambiguities in the choice of observables and hopping rates, since the natural observable is the 
position : we are interested into the response of the position to an external bias, and in the thermal fluctuations of 
the position. 



II. STOCHASTIC DYNAMICS OF THE TRAP MODEL 



In this Section, we present the Master Equations that describe the stochastic dynamics of the trap model in the 
absence and in the presence of an external field, and recall some useful notions and properties. 



A. Trap model without external field 



1. Master Equation in a given sample 



The dynamics of the one-dimensional trap model is defined by the following Master Equation for the probability 
Pt(i) to be at time t on site n |26j 

dP t (n) _ P t (n + l) | P t {n-l) P t (n) 
dt 2r n +i 2t„_i r n 

supplemented by the simple initial condition at t = 

Pt=o(n) = S n>0 (2) 

This means that after arriving on site n, the particle spends on site n a random time t distributed with the exponential 
distribution of mean r„ 

t_ 

U n {t) = -e~Tn (3) 

and then jumps with equal probability (1/2) to one of its two nearest-neighbor sites (n — 1) and (n + 1). 

The trapping times {t„} are quenched random variables in the Master Equation and represent the disorder : at 
each visit to the site n, the particle is submitted to the the same mean trapping time r„. A given sample corresponds 
to a particular realization of the trapping times on the whole line {..., T_fc, ...,tq,ti, ...,7%, ...}. 



2. Law for the disorder 



The trapping times are usually taken to be independent random variables distributed with a law presenting the 
algebraic decay jU 

q(r) oc (4) 
This choice comes from the interpretation of the trapping times as Arrhenius times associated to energies E n 

r n = (5) 
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where the distribution of the energies E n presents the exponential tail 

p(E) cx e - ^ (6) 

to mimic the lowest energies statistics found in the Random Energy Model |37j and the distribution of free energy 
of states in the replica theory of spin-glasses E3 ■ This choice is moreover reinforced by the exponential tail of 
the Gumbel distribution which represents one universality classes of extreme-value statistics. In this interpretation, 
the exponent /i in J3J is directly related to the temperature /i = T/T g and T g represents a glass transition. Indeed 
for /it < 1, the mean trapping time J d,TTq(r) is infinite and this directly leads to aging effects. In particular, the 
characteristic length scale £ (t) for the position n at large time t follows the sub-diffusive behavior obtained by various 
methods 



■istic length scale £( 
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m ~ 1— (7) 

for a < 1 and not the usual scaling in t 1 / 2 of the pure random walk. For more detailed results concerning the 
localization properties, the averaged diffusion front as well as aging and sub-aging properties, we refer the reader to 
the recent works [3l|, |32| ■ 

B. Trap model in the presence of an external field 

In the presence of an external force field /, the master equation becomes E3 

^-M = P (f) (n + 1} |_JL + P (/) (n _ l)| — - - P^in f- 3 +6 (8) 
at £T n+ i 4T n -i iT n 

where the hopping rates have been chosen to satisfy the detailed balance condition 

e -0U(n) W U) -PU n+lv Af) (Q) 

{n— m-fl} {n+1— mj \ ' 

where the total energy U n contains both the random energy (—E n ) (0 of the trap n and the potential energy (— fn) 
linear in the position n induced by the external applied field / 

U n = -E n - fn (10) 

The condition guarantees that in finite sizes, the dynamics will converge towards the steady state given by 
Boltzmann equilibrium P eq (n) ~ e~ P u \ n > . 

III. CONSERVATION OF THE MEAN POSITION IN ANY GIVEN SAMPLE 

As already emphasized in |33j |. the trap model has a very special property : when there is no external field / = 
the mean position defined by 

<n £ nP t U=0 \n) (11) 

71— — OO 

is a conserved quantity in any given sample ! Indeed, the master equation directly yields 

d < " {t) = (12) 

for any realization of the trapping times {r„}. This comes from the fact that the hopping rates only depend on the 
initial state, and not on the final state. As a consequence, when a particle gets out of any trap, it always jump 
with equal probability (1/2) on the right or on the left. So the trap model possesses a local symmetry for the jump 
probabilities in the two directions. This is in contrast with other models of random walks in random media, such 
as for instance the Sinai model, where this local symmetr y d oes not hold, and where as a consequence, there is a 
non-trivial mean position in each sample < n > f—o (t) ^ |41| , and it is only after the average over the samples that 
the symmetry in the two directions is restored, i.e. < n >/=o (t) = 0. 

The property (|12f> is all the more remarkable that the disorder completely breaks the symmetry n — > — n in any 
given sample! In particular, the diffusion front Pt(n) in a given sample is of course not symmetric in n — > — n, and it 
is only the average over the disordered samples that is symmetric in n — > —n, because to each realization {t„}, one 
may always associate the reversed realization {t_„}. 
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IV. NON-LINEAR RESPONSE TO A BIAS APPLIED FOR t > 
A. A non-linear Fluctuation Theorem 



We now consider a given sample of the trap model and we wish to compare the dynamics in the presence of an 
external field (— /) with respect to dynamics in the presence of an external field (+/). It is clear from the two 
associated master equations iJSJ that we have the remarkable simple property 



>(+/) 



(n) 
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We may thus express the differences of the mean position for the two cases as 



+00 



+ OC 



< n >+/(*)-< n >_/(*) = n [ p t +f \ n )- p t~ f \n)] = ^ n{\ - e~^ n )P^ f) {n) 



(13) 



(14) 



B. Fluctuation-Dissipation Relation in the linear response regime 



A first consequence of ((14(1 is that the susceptibility with respect to the force / reads 



d <n> f (t) 



df 



I < n 2 > /=0 (t) 



(15) 



and thus the Fluctuation-Dissipation Relation is valid for fixed t in the limit / — > in the form of the Einstein relation 



< n >/ (t) £ o ^ < n 2 > /=0 (t) 



(16) 



Breakdown at long times of the linear response 



A second consequence of 1(14(1 is that for fixed /, the linear regime 116(1 of the response is valid only if the linearization 
of the factor (1 — e~@* n ) is consistent, and thus it yields the following criterium for the validity for the linear response 



/3/£ (/ =o)W«l 



(17) 



where &f—Q)(t) represents the length scale for the displacement at time t in the unbiased trap model and follows the 
anomalous scaling JJJ. The linear response is thus valid only up to a characteristic time presenting the following 
dependence in the force / 



(18) 



D. Consequence for the disorder averaged diffusion front 

A third consequence of 1(14(1 is the following explicit form for the symmetry breaking between (n) > and (— n) < 
for the disorder averaged diffusion front in the presence of the external field / > 



py\-n)=e- 0fn P t if) (n) 



(19) 



To obtain this property, it is convenient to associate to each realization T> = {r„} of the disorder the reversed 
realization T> R = {t^} simply defined by t„ = t_„ for any n. It is then clear that the the dynamics in the presence 
of the external force (+/) in the sample T> in the coordinate (n) is completely equivalent to the dynamics in the 
presence of the external force (— /) in the sample T> R in the coordinate (—n), i.e. more explicitly we have for the 
diffusion fronts 



pV\n,t)=p(- R f) (-n,t) 



(20) 
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So the relation l|13[) between the diffusion fronts in the same disordered sample in the presence of two opposite external 
forces (+/) and (— /), may also be rewritten as a relation between the diffusion fronts in the two reversed disordered 
samples in the presence of the same external force (+/) 



4 +/ W) 



(21) 



and thus the average over the samples D immediately yields the simple property (|19(1 . As a comparison with other 
models of random walks in random media, we may cite the Sinai model : the explicit expressions given in |4lj for the 
diffusion front in the presence of a bias show that the property (|19|l is not satisfied for the Sinai model. 

V. NON-LINEAR RESPONSE TO A BIAS APPLIED FOR t > t w 

We now consider the following aging experiment |33| , which is very usual in the study of disordered systems 0, El 
: the system first evolves with no external force f — during the interval [0,^], and then evolves in the presence of 
the external force (/) for t > t w . We consider the two-time diffusion front P^' tw '(n, t; n w , t w \0, 0) which represents 
the joint probability to be at position n w at time t w and at position n at time t (where it is assumed that t > t w ). 
Since the dynamics takes place with no bias up to t w , it is clear that for any sample, we have from the property 1|12|) 

<n w > Uw =YjY^ n W P U '' tw) (n, t; n w , t w |0, 0) = (22) 

n n w 

so that the mean position at time t also corresponds to the mean displacement between t w and t 

<n>f, tw (t) =Y,J2 nP(f '' ta) ( n > t > nw ' tw \°'ty =< ( n-n ™) >/.*« (*) ( 23 ) 

n n w 

Again, we wish to compare the dynamics in a given sample for the two cases (+/) and (— /). The generalization of 
(H2J reads 

p(+/;M( M;7Wtu | 0i 0) = f[n _ nw) 

P(-f^)(n,t;n w ,t w \0,0) 1 ' 

which leads to the following property for the relative mean displacement between the two times t w and t 

<(n-n w )> +f]tw (t)-<(n-n w )>-f., tw (t)= £ £ (n ~ n w )(l ~ e-^ n ' n ^)P^ f ^\i,t; lw ,t w \0,0p5) 



n— — oo n w — — oo 



A. Fluctuation-Dissipation Relation in the linear response regime 

So for fixed times (t,t w ), there is a linear response at small field / given by 

Bf 

< (n - n w ) > f , ta (t) £ q ^<(n- n w ) 2 > /=0 (i) (26) 
which generalizes the one-time Einstein relation l|16() . 

B. Breakdown of the linear response regime 

But for fixed waiting time t w and fixed external field /, the validity of the linear response is limited in time t by 
the condition of the validity of the linearization of the factor (1 — e"^™ - ™ 1 ")) in (|25|l which reads 

Pftf=o{t,t w ) « 1 (27) 
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where £/ = o(i, t w ) represents the length scale for the relative displacement (n — n w ) between the times t w and t in the 
unbiased trap model. In the aging regime where t w oo, t — > oo with the ratio (t — t w )/t w > 1 is fixed, the relative 
displacement scales as [^3, 

C/=o (*»*») ~ (*-*«;) ^ (28) 
The validity of the linear response regime is thus again limited in the aging regime by 

t-t w «t^{f) (29) 

in terms of the characteristic time Ijl8|l . 

C. Consequence for the two-time disorder averaged diffusion front 

The generalization of the relation (|2(J|) between the diffusion fronts in two reversed disorder samples T> and T> R to 
the two-time aging experiment reads 

P^ f;tw \n, t; n w ,t w \0, 0) = P^ /; * lu) (-n, t; -n w , t w \0, 0) (30) 

So the relation (|24JI between the diffusion fronts in the same disordered sample in the presence of opposite external 
forces (+/) and (— /), may also be rewritten as a relation between the diffusion fronts in the two reversed disordered 
samples in the presence of the same external forces (+/) 

Py' w \ n it'i n w ,t w \0, 0) _ e pf( n - nw ) 
P v f A tw) (-n,t;-n w ,t w \0,0) 

The average over the samples T> thus yields the following generalization of (|I9fl for the asymmetry in n induced by 
the external force / on the two-time disorder averaged diffusion front 



p(/)(_ n , f, -n w ,t w \0, 0) = e-W n - n ^P(f)(n, t; n w , t w \0, 0) (32) 

VI. GENERAL FORM OF THE SPECIAL SYMMETRY LEADING TO A SIMPLE NON-LINEAR 

FLUCTUATION THEOREM 

In this Section, we investigate the precise condition under which a general master equation will yield a simple 
non-linear Fluctuation Theorem such as the one found for the trap model. 

A. Master Equation with detailed balance in the presence of an external field 

We consider a general Master Equation, where the configuration of the system is denoted by C 

dP ^ C) = ^ P t if) (C')W^ (C -h. C) - P t if] (C) Y, W {f) (C - C) (33) 

where the hopping rates in the presence of an external field / linearly coupled to the observable b are chosen as 

WW(C -» C) = W°(C -> C r )e^ b ^- b ^ (34) 
to satisfy the detailed balance condition 

W U)(C'^C) W°(C'^C) K ' 

We consider the simple initial condition 

P t ( ll(C) = 5c,c (36) 

where (Co) is some arbitrary given configuration. 
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B. Symmetry condition for the total leaving rates 

Let us now define the auxiliary function 

Q t {C)=Py\C)e^^- h ^ (37) 

It satisfies the initial condition 

Qt=o(C) = S CtCo (38) 

and evolves with the master equation 

- £ Qt(C')W (C - C)e^ b ^- b ^ - Q f (C) £ W°(C - C')e¥[*(c')-KC)] (39) 
c c 

We wish to compare with the master equation in the presence of the external field (— /) which reads 

( — /) 

^ m =J2 P t~ f \ C ') W °( C ' "> C)e^ b ^- b ^ - Pt f) (C)J2W {0 HC -» C ')e^ [fc(C) - fc(C,)1 (40) 

* a c 

The "arriving" hopping rates in the configuration C are thus identical in (|39() and l|40l) . but the "leaving" rates out 
of the configuration C are obviously different term by term. The important quantity is actually the total leaving rate 
out of a configuration C in the presence of the external field (+/) which reads 



w 



(+/) (C) = W°(C -» C^WO-HO] (41) 



07/ / 



which is of course a priori different from w ou l (C). This means that generically, the local energy landscape around 

the configuration C will be distorted in different ways by the application of (+/) and by the application of (— /). 

( — f ) 

So the condition to have the simple relation P t — Qt, leading to the following non- linear Fluctuation Theorem 

P t U) ( C ) _ J3f[b(C)-b(Co)] U9) 

is the presence of the special symmetry 

w£P \C) = w^P \C) for any configuration C (43) 

This symmetry is satisfied in the trap model for any configuration in spite of the presence of the disorder, because 
the hopping rates of the unbiased case depend only on the initial configuration and not on the final configuration 

W°(C -> C) = W°(C) (44) 

and thus the disorder contained in the hopping rates disappear from the symmetry condition 143|) . which is then 
satisfied because the observable b which is the position can increase or decrease by one with equal probabilities. 

C. Symmetry condition in the linear regime 

At linear order in the external field /, the symmetry condition l|43|l simplifies into the following condition for the 
unbiased hopping rates 

Yy>(C')-b{C)]W°{C ->C) =0 for any configuration C (45) 
c> 

This exactly means that the mean value < b > is a conserved quantity for the unbiased dynamics and for an arbitrary 
initial configuration 



d < b > 

dt ~ dt 



^y>{C) ~b{C)]W^\C ->C) 



(46) 



8 



which vanishes for an arbitrary set of probabilities P} (C) only if the condition (|45|l is fullfilled. For the trap model, 
this condition corresponds to the conservation of the mean position in any given sample l|12f) . 

If the condition 145|) is satisfied, then the Fluctuation-Dissipation Relation is valid for fixed t in the limit / — > in 
the form of the Einstein relation 

< (6 - b ) > f (t) £ Q & -l < (b - 6 ) 2 >/=o (*) (47) 



D. Discussion 

It would be interesting to find other examples of stochastic dynamics, which satisfy the dynamical symmetry l|43|l . 
apart from the trap model considered here. A simple way to satisfy the condition (|43[) would be to be able to associate 
to each transition C — > C an other transition C — » C" which has the same transition rate in zero-field W^°\C — > C") = 
VF (0) (C -> C) and such that the variation of the observable b is exactly the opposite b(C") ~ 6(C) = -(6(C) - 6(C)). 
It is however a very strong property, since it should be true for any starting configuration C, and we have not been 
able to find another example in the most standard dynamics usually considered in statistical physics. 

VII. CONCLUSION 

In this paper, we have obtained that the stochastic dynamics of trap models has a very special symmetry that 
implies a simple non-linear Fluctuation Theorem for the one-time and two-time diffusion fronts in any given sample. 
We have obtained the following consequences : (i) in the linear response regime, the Fluctuation-Dissipation Relation 
or Einstein relation (|26|l is valid even in the aging time sector; (ii) for a fixed waiting time t w and fixed external field 
/, the validity of the linear response regime is limited in the time (t — t w ) by the characteristic time t^{f) ~ (/3/) _T +^; 
(iii) in the non-linear regime, we have obtained very simple relations for the asymmetry in the position n induced by 
an external force on the one-time and the two-time disorder averaged diffusion fronts. 

The results (i) and (ii) are in full agreement with the scaling analysis and the numerical simulations presented in 
|33| . and we hope that our approach in some sense "explains" why aging trap models, which are known to present 
str ong localization effects and to be completely out-of-equilibrium (see for instance the effective dynamics described 
in |32|h nevertheless satisfies the Fluctuation-Dissipation Relation in the whole linear response regime, in contrast 
with many disordered systems which present a non-trivial violation in the linear response regime |& |g, |2j - 

In a future work, we will consider again the same problem of the effect of an external force on the trap model with a 
complementary point of view : here we have used a special symmetry of the Master Equation to derive general results 
which are thus valid in the whole aging phase < /i < 1 J3J ; in the forthcoming paper j^, we will obtain various 
exact explicit results characterizing response properties in the high disorder limit fi — > 0, by generalizing the previous 
real-space renormalization approach for the unbiased trap model |32( to include the presence of an external bias. In 
particular, we will show that the explicit expressions for the averaged diffusion fronts indeed satisfy the constraints 
(iii) on the asymmetry in n obtained here by a general argument. 
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